Hadron Spectrum and Decay Constant from Np = 2 Domain Wall QCD 

Taku Izubuchi a < b [RBC Collaboration] * 

a RIKEN-BNL Research Center, Brookhaven National Laboratory, Upton, New York 11973, USA 
b Institute for Theoretical Physics, Kanazawa University, Kanazawa 920-1192, Japan 

We report on the first large-scale study of two flavor QCD with domain wall fermions (DWF). Simulation has 
been carried out at three dynamical quark mass values about 1/2, 3/4, and 1 m s t ra nge on 16 s x 32 volume with 
L s = 12 and a" 1 « 1.7 GeV. After discussing the details of the simulation, we report on the light hadron spectrum 
and decay constants. 



1. Introduction 



jectory, scaled by V and the MD step size, At, 



We have continued to investigate two flavor 
QCD using dynamical DWF. We present results 
of the hadron spectrum and the decay constants 
as a project of the RBC collaboration. For other 
quantities on the same dynamical ensembles, sec 
other RBC contributions in these proceedings. 

In this work we choose the DBW2 gauge action 
at f3 = 0.80, in which the negative coefficient to 
the rectangular plaquette suppresses dislocations 
of configuration and drastically reduces the resid- 
ual mass, m res , for relatively small L s ~ O(10) 
in quenched simulation 1 . We simulated three 
sea quark mass m sea = 0.02,0.03, and 0.04 using 
Np = 2 domain wall fermion with fifth dimen- 
sional extent L s = 12 and domain wall height 
M 5 = 1.8. L s = 12 



2. Ensemble generation 

Tabled summarizes the HMC-$ evolution. An 
interesting result is that the acceptance, P acc , 
is insensitive to m sea . More precisely, at least 
for our particular run parameters: m sea > 
fn strange / ^2, the squared energy difference be- 
tween the first and the last configuration in a tra- 
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stays same for all sea quark masses in the simula- 
tion, while an empirical estimation for staggered 



or Wilson fermions would be Cah 



,a~2, 



in which case number of steps has to be in- 
creased proportional to mj^ a to keep the accep- 
tance P acc w erfc(i/ ((AH) 2 ) /8) constant. By 
introducing the new force term described in [2], 
Cah is reduced to ~ 40 % of the old one. 

We have implemented the chronological in- 
verter technique of 3 , in which the starting vec- 
tor of the conjugate gradient algorithm (CG) in 
a MD step is forecasted by the solutions in the 
previous steps. Nqq in the table is the number 
of matrix multiplication in CG without forecast- 

(7) 

ing, while Nqq is the count using the previous 
seven solutions. The number seven is about the 
point where the precision of the forecast is satu- 



rated. N, 



(tot) 
CG 



is the total number of CG count 



in a trajectory. By a scaling ansatz, N^ G = 
C'i(m sea + m res )~^, fa ~ 1 and (3 7 — 1.5 from 
the table. Although fa > fa, the coefficient, d, 
is much smaller for forecasting and thus N^q^ is 
reduced roughly by factors of two to three in the 
simulation points compared to the case without 
the forecast. 

Ti n t in the table is an estimation for the inte- 
grated auto-correlation length for 1 x 1 plaquette 
using whole trajectories by the truncated sum at 
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Table 1 

The details of L s = 12, M 5 = 1.8 DWF HMC-$ evolution on 16 3 x 32 lattice with (3 = 0.80 DBW2 
action. 
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50 with an error from jackknife blocks of length 
100. Tint are all equal within the quoted errors. 
The autocorrelation length of axial current corre- 
lator measured every ten trajectory is < 50 tra- 
jectories 2 

Although these observations are encouraging 
for DWF simulations with lighter quark masses 
in the future, we note that they may only be true 
for relatively heavy m sea . 

3. Physical Results 

For each of three 

wisea i we measure on 94 lat- 
tices separated by 50 trajectories, leaving out the 
first ~ 600 trajectories to allow the evolution to 
thcrmalized. All quoted errors are from jackknife 
estimate of the statistical error, correlated fit to 
hadron propagators is made using a single covari- 
ance matrix computed on the entire ensemble. 

By a linear diagonal extrapolation, m sea — 
rrivai — ► 0, we found m res = 0.00137(4) or < 5 
MeV, which is larger than quenched DBW2 val- 
ues for the same L s and a , but it is still an order 
of magnitude smaller than the input quark mass. 
The impacts of m res on the operator mixing for 
Bk is discussed elsewhere |4I5) . 

From non-gauge-fixed wall-point pseudo-scalar 
correlator, (J5 Jg ) we extract pseudo-scalar decay 
constant: 



f P s = — o-(m/ + m res ) (0| J5 \w,p = 0) 



in 



(2) 



and linearly extrapolate to the chiral limit, 



m lal = m lea 

m res , m q sea — m va i + m res , to obtain an esti- 
mation for its chiral value, / = 0.078(1). The 
fit to the next-to-leading order (NLO) partially 



with 



2 The topological charge and the chiral condensation (with 
many random hits) show longer autocorrelation times. 



quenched chiral perturbation theory (PQChPT) 
formula [6l7| for f ps did not describe the data 
in the neighborhood of the fit range too wellp]. 
However, the following results depend on the 
value of / only mildly, so we will use the esti- 
mation from the linear fit. 

Fitting Coulomb gauge-fixed wall-point corre- 
lators, we obtain the mass of the pseudo-scalar 
mesons and the vector mesons. It is worth men- 
tioning that a simple linear extrapolation of 
t° m lai sea = is zero within the statistical er- 
ror. This did not happen in the quenched case: 
to 2 s — at to/ « — (2 — 3) x m res pQ. This 
is consistent with the difference between chiral 
logarithms (rn 2 vs jm q ~ 2Bo + cm q log to 9 ) vs the 
quenched one [rr? ps jm q ~ logm 9 ) at small to 9 . 

The NLO PQChPT formula, 
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is used to fit the pseudo-scalar mass obtained 
from axial current correlators and pseudo-scalar 
correlators for various mass ranges. The NLO fit 



is constrained as to 2 
i 



at m va i = — TO r 



as it 



must due to the universality of the low energy do- 
main wall fermion theory, and / estimated above 
is used as an input. The results of fit of the fit are 
shown in Figure and in Table [21 Note that Bg 
is stable under change of fitting range, while the 



3 



Li coefficients are not. When the heavier mass is 
included in the fit, the \ 2 /dof increases. 

Table 2 

Parameters from chiral perturbation theory fits 
to the values of at to/ = m sea . va i < nrSp 11 ^ 
computed from the pseudo-scalar wall-point (up- 
per two column), and axial- vector wall point. \ 2 
is from uncorrelated fits in m/. Li refer to Gasser- 
Leutwyler low energy constants multiplied by 10 4 
at A x = 1 GeV. 

mf ax) rYdof 2 go L 5 -2L 8 L i -2L e> 

0.03 0.1(1) 4.0(3) -1.5(7) -2(1) 

0.04 2(1) 4.2(1) -0.2(4) -1-1(4) 

0.03 0.3(2) 4.0(3) -1.9(8) -1(1) 

0.04 1.9(9) 4.2(1) -0.4(4) -0.8(3) 



mula for non-degenerate quarks are determined 
by the fit of the degenerate quark meson. It is im- 
portant to note that to and m s as defined above 
are bare quark masses; the renormalized quark 
mass is defined as Z m m q = Z m (m + m res ), where 
Z m is a scheme and scale dependent renormaliza- 
tion factor [3]. 

Extrapolating the measured decay constant at 
the partially quenched points (m va ^ sea < 0.04) 
using the linear ansatz to to and m s , we find f n = 
134(4) MeV, f K = 157(4) MeV, and f K /U = 
1.18(1), which agree better with experiment than 
quenched DWF simulations. 

Sommer scale, r$ (with ro = 0.5fm ), from 
the static quark potential gives almost the same 
value for the lattice spacing 8 . A similar analy- 
sis has been carried out for baryons, and the chi- 
ral limit of the linear extrapolation gives Mn/M p 
= 1.34(4), which is larger than the experimental 
value. This would be expected from the relatively 
small spatial volume. 
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Figure 1. vn^J (m va i+m res ) as a function of m va i 
for each m sea = 0.04,0.03,0.02. The curves are 
from the results of NLO fit using m sea ^ va i < 0.03. 

A linear extrapolation of the three vec- 
tor meson mass points m va i = m sea , yields 
a -1 =1.69(5) GeV at the quark mass, to/ = 
to = 0.0002(1) corresponding to the neutral 
pion mass using ©. Similarly, but using the 
NLO ChPT prediction for non-degenerate valence 
quarks, m s = 0.045(3) is obtained from the phys- 
ical kaon mass. The coefficients in the NLO for- 



4. Conclusion 

We have generated configurations of lattice 
QCD with two-flavor dynamical DWF at a -1 « 
1.7 GeV, with a spacial volume V « (1.9fm) 3 . 
From the ~ 3 x 5000 trajectories in this work 
obtained with light dynamical quarks and small 
residual quark mass, TO stra „ ge /2 m sea + m res 
^strange or 0.54 ^ m ps / 'my ^ 0.65. We have 
tried the NLO ChPT fit for the mass and a sim- 
ple linear fit for the decay constant. The results 
for the physical decay constants are closer than 
those in quenched DWF simulation. For further 
details of the simulations and results, we refer the 
reader to the forthcoming paper 5 . 
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